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Fermionic superfluids can undergo phase transitions into different kinds of normal regimes, loosely 
characterized by whether Cooper pairs remain locally stable. If the normal phase retains strong 
pairing fluctuations, it behaves like a liquid of vortices, which has been observed in cuprate super- 
conductors. We argue that analogous strongly correlated normal states exist in two-dimensional 
neutral fermion liquids near unitarity, where superfluid is destroyed by fast rotation. These states 
have non-universal properties, and if they develop as distinct thermodynamic phases they can be 
characterized as quantum Hall states of Cooper pairs. The formal analysis is based on a model 
with SP(2A'') symmetry that describes the quantum critical region in the vicinity of a broad Fesh- 
bach resonance. We explore the pairing phase diagram and demonstrate that the considered model 
has macroscopically degenerate bosonic modes in the normal phase, to all orders in 1/A'^. It takes 
finite-range interactions to lift this degeneracy, making the Abrikosov flux lattice of the superfluid 
particularly susceptible to quantum melting. 



Recent years have seen many developments regarding 
strongly correlated normal phases of fermionic supercon- 
ductors and superfluids. It is becoming increasingly clear 
that the unconventional normal phase of cuprate super- 
conductors results from quantum fluctuations that de- 
stroy the long-range phase coherence of Cooper pairs 
These fluctuations can be ascribed to the motion 
of vortices, and indeed the normal state has properties 
of a vortex liquid [1, 01 • A byproduct of the fluctuations 
and crystal lattice is the emergence of density modula- 
tions, [4 j_ al, which are nowdays observed in great detail 

Ultra-cold atoms, being much simpler than electronic 
systems like cuprates, are very attractive for studies of 
many-body physics This paper explores the possi- 
bility of obtaining analogous strongly correlated normal 
states using neutral cold atoms in continuum. We con- 
sider two-dimensional two-component fermion gasses in 
which the superfluid phase is destabilized by fast rota- 
tion at low temperatures 12, 13, 14 1. We argue that the 



unconventional normal states are most likely to be stable 
near a broad Fcshbach resonance, when two-body scatter- 
ing length is very large. Due to the absence of long-range 
forces between neutral particles, a large scattering length 
is what provides interactions potent enough to stabilize 
strongly correlated normal states. Experimental progress 
toward achieving such states is still limitted, but the over- 
all pace of development in the field gives many reasons 
for optimism 1^, 3, 17|. 

Static three-dimensional cold-atom systems near a 
broad Feshbach resonance have striking universal prop- 
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erties 
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20|. This universality is controlled by 
a quantum critical point at the zero density Feshbach 
resonance flit , which defines the unitarity limit. The su- 
perfluid in this regime interpolates between the Bardeen- 
Cooper-Schrieffer (BCS) state of Cooper pairs and Bose- 
Einstein condensate (BEC) of diatomic molecules, and 



insulating states in optical lattices interpolate between 
the two analogous limits of band and Mott insulators re- 
spectively [2ll |. 

The situation is different in two-dimensional rotating 
cold atom systems because some universal properties are 
lost. Nevertheless, the unitarity limit is well defined be- 
cause it originates in the nature of interaction potential 
between atoms. Consider a theory of 2N fermion species 
(a € {Ti}i * = I-- - ^) coupled to a bosonic Cooper pair 
field $ in d dimensions {h = 1): 



S = 



(1) 



This is a generalization of the "two-channel" model, 
which we recover by setting iV = 1. The operator H^'^' is 
proportional to detuning from the Feshbach resonance 
(a part of Eq. ([9])). In the absence of rotation at T = 0, 
/i = (zero density) and = (Feshbach resonance) this 
theory is at a fixed point of renormalization group (RG) , 
and no o per ators other than the written ones are relevant 
in d > 2 [l9|. Therefore, only the zero-range interactions 
matter, and finite short-range forces have negligible im- 
pact on macroscopic properties near the fixed point. Note 
that no particular scale characterizes interactions, so that 
the usual perturbation theory in interaction strength is 
not applicable. For that reason we generalize to large 
N and obtain controlled approximations by perturbative 
expansions in powers of 1/A^. 

We now introduce rotation at angular velocity co and a 
potential well to make the system quasi two-dimensional: 
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where V x A = zB = zmujc = z2mu}. This action 
describes dynamics in the rotating frame, assuming that 
centrifugal forces are exactly balanced by a trap. We 
choose this simplification in order to focus on many-body 
effects without being side-tracked by trap effects. 

The main utility of considering the model ^ comes 
from its simplicity and striking properties that character- 
ize the unitarity limit. We shall reveal an unusual pairing 
instability in this model, which involves a macroscopic 
number of soft modes that are degenerate to all orders in 
An implication is that all possible arrangements of 
vortices in the superfluid differ very little by free energy 
close to the pairing transition, so that small perturba- 
tions can easily mix them and melt the vortex lattice. 
Note that the phases of the order parameter are only al- 
gebraically correlated in the two-dimensional superfluid 



with a vortex lattice [H,!!!,!!!. 

The model ([2]) contains only zero-range interactions 
(the Hubbard-Stratonovich field <E> has no bare disper- 
sion at the fixed point of ([T])). We demonstrate using 
RG that finite-range interactions between fermions are 
perturbations to ^ that can shape phases and induce 
transitions. Physically, this comes from the fact that the 
bare fermion spectrum in ([2]) consists of dispersionless 
macroscopically degenerate Landau levels with energies 
e„ = niOc — (we measure energy with respect to the 
redefined chemical potential fj,' = ii — Wc/2). The lack of 
dispersion reduces dimensionality in RG and makes even 
arbitrarily weak perturbations extremely potent in lifting 
the degeneracy. We set up RG by generalizing ^ to d 
dimensions in the Landau level basis using the Landau 
gauge, and include additional allowed terms: 
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We suppress spin a and fiavor i indices for brevity. Since 
the bare fermion states are localized in the plane per- 
pendicular to the axis of rotation, it is appropriate to 
not rescale kx in RG. At zero temperature and density, 
fermion self-energy vanishes so that all renormalization 
comes from the boson field $. The RG fiow equations 
can be calculated exactly to all orders of perturbation 
theory when U2 and the remaining omitted couplings are 
zero, since then the renormalization of $ involves sum- 
mation of a geometric series of bare fermion bubble dia- 
grams. Following the procedure in Ref . [l9| . we find that 
under rescaling r'^ = e~^r±, t' = e~^'T, tp' = e^''/^~"'^-''-0, 
$' = e^''/^^^^'^ the exact fiow equations are: 

If-"' 

where a and b are cut-off dependent constants whose val- 
ues are not important for the present discussion, and 

n(o) 

cx ly. The interacting fixed point at ji' ~ ujc = i' ~ 0, 
g* = 1/(3 — d/2)/{hN) defines the unitarity limit in any 
dimension d, and shows how perturbation theory can be 
justified for large N . Appart from /x', ujc and various 
couplings u„, which may include non-local short-range 



potentials and hence multi-particle collision terms, can 
be relevant at this fixed point. For n-particle scattering 
Un we find sufficient indication at the tree-level: 



-^ = [d+{2-d)n]ur, + 0{u') (5) 

irrespective of the spatial dependence of interaction po- 
tentials in the plane perpendicular to rotation axis, since 
coordinates do not rescale in this plane. Therefore, all 
Un are relevant in d < 2, while in d > 2 they are relevant 
for n < d/{d — 2), with an upper bound n N beyond 
the tree- level for > 1 and d>2. 

The physical consequences of this scaling will be ap- 
parent shortly. However, we first elucidate the effects 
associated with unitarity by ignoring the perturbations 
{un 0) and absorbing g in the definition of Hence, 
we now focus on the simple theory ([2]) in two dimensions. 

Let us start in a quantum Hall state of unpaired 
fermions and consider the pairing instability. This in- 
stability shows in the poles of the boson field Green's 
function given by the inverse fermion bubble diagram in 
particle-particle channel, G$ = (A^H)"^. The bubble di- 
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agram in Landau representation is (without n*^"^): 
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Here, as in ([3]), the Landau gauge quantum numbers n 
and /ct determine fermionic and bosonic wavefunctions: 
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where Hn{x) are Hermite polynomials. Note that boson 
"charge" is twice that of fermions. The vertex function 
in Landau representation is readily derived from ([7]): 
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Note that the vertex does not depend on the boson mo- 
mentum px, as a curious consequence of the Landau level 
degeneracy. The bare fermion Green's functions in Lan- 
dau representation is simply Gn{iuj) — {—iuj + en)~^- 

Pairing instability can be directly found from the bub- 
ble diagram matrix H at zero Matsubara frequency as the 
onset of a first negative eigenvalue. The matrix elements 
of n in Landau representation are given by Q . We note 
that ([6]) is not gauge invariant, but the spectrum of H 
is: gauge transformations are manifestly unitary trans- 
formations of n in position representation. Therefore, 
we can safely fix gauge and look for the instability in any 
convenient representation. 

Before preceding with calculations, it is necessary to 
regularize the ultra-violet divergent expression ([6]). In 
the absence of rotation, regularization leaves behind a 
chemical potential term for the $ field, proportional to 
the detuning v from the Feshbach resonance The 
same procedure is problematic in the rotating case be- 
cause Landau orbitals are localized and scattering length 
a = whose formal definition rests on the scatter- 

ing of free particles, is not entirely meaningful. Only for 
a Rc, the cyclotron radius at Fermi energy, v retains 
its physical meaning. We formally extend the definition 
of detuning from this limit to all values v. Using the re- 
lationship between a and two-body scattering vertex, we 
regularize U by adding "AH = 0" to ([6]): 
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where is a quantum well confinement length scale in z- 
direction and n(0) = (p^ = 0|n|p^ = 0) is the projection 
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FIG. 1: (color online) (a) Critical temperature as a func- 
tion of chemical potential for several values of detuning 1/ = 
1/ / (imcoaz): the outermost curve is BEC limit u = —0.9, the 
innermost is BCS limit u — +0.9, unitarity v — is dashed 
(AJ- = 0.1). The vertically stretching lines are normal state 
constant density contours, with increment 0.25B/{2tt). (b) 
through (d): The evolution of the transition line with detun- 
ing on the BCS side {u is 0.4 in (b), 0.5 in (c), 0.6 in (d). 



of n to plane- waves at zero momentum: 

nrm-Vn i^n ~ i)ii{2n' ~ ly.i 

The expression ^ guaranties that the ultra-violet di- 
vergence is cancelled in all matrix elements, without in- 
troducing arbitrary (unknown) features at cut-off scales. 

Figure [1] shows pairing second-order phase transitions 
in the limit N ^ oo. The deeper one goes into the 
BEC regime, the faster the growth of critical temperature 
with chemical potential (density). In the BCS regime, 
however, the paired phase breaks up into dome-shaped 
islands at low densities, obtained when chemical poten- 
tial crosses a fermionic Landau level. Similar phenomena 
have been proposed in cuprates [25[ . Once a paired state 
becomes a connected region at sufficiently large densi- 
ties, there are trailing islands of unpaired states sitting 
between the Landau levels. Therefore, unpaired normal 
states can exist at T = 0. 

We now establish a crucial fact: ^ is independent 
of Px to all orders of 1/iV. First, the written lowest or- 
der term in the expansion involves fermion energies 
that do not depend on ±kx + Px/2, and vertices ([8]) 
that do not depend on px- This term defines the bare 
boson propagator G$ cx iV~^, using which we generate 
higher order corrections of ([5]). Let us apply the follow- 
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FIG. 2: Feynman diagram labeling examples. 



ing labeling rules: each fcrmion propagator shall carry 
momentum kx +Px/'2 in the arrow direction, while each 
boson propagator shall carry Qx + Px (see Fig 121). The 
transfer of px is automatically conserved. Each added 
vertex takes momenta kxi + Px/'^ and kx2 + Px/'^ at its 
fermionic terminals, but there is no dependence on px 
since only the difference of momenta at the fermion ter- 
minals matters in ([5]). The added bare boson lines take 
Px contributions, but according to ^ they do not de- 
pend on the transferred momentum. Therefore, order by 
order, no dependence onpx is introduced in any Feynman 
diagram. 

This means that the poles of the exact boson Green's 
function are macroscopically degenerate in the normal 
phase, having no dependence on px- Pairing instability 
involves an infinite number of bosonic modes going soft 
at the same time. What are the physical consequences? 

Paired states are described by order parameters $(r) 
corresponding to static arrangements of vortices whose 
total number is AB/tt, A being the system area. In 
general, free energy density J- is minimized by an order 
parameter with a periodic array of vortices: 

$(r) = B-iJ2J2 J2 ^nMi+nySqj{r) , (11) 
n ;=o j 

where the integer Uy and momentum scale dq de- 
termine the flux lattice periods Aa; — 2'K/6q and 
Ay = ny6q{2B)~^ (note that <i>(r) has aperiodic gauge- 
dependent phase). The amplitudes 4>n,i, as well as Sq and 
Uy, are obtained by minimizing JF, and in normal circum- 
stances correspond to the triangular Abrikosov lattice: 
Sq — (27r\/3i?)^/^, Uy — 2, (f>o_i — i^o.o- The minimum 
of is well defined and unique deep inside the paired 
phase. However, near the pairing transition a large num- 
ber of different paired states become rapidly competitive 
by free energy because: 1) the evolution of with action 
parameters is smooth and 2) J^[$(r)] approaches J^[0] 
at the same rate for all $(r) owing to the degeneracy of 
n = d^JF/d<i>(r)^. As an illustration we compare the free 
energy minimums of two flux lattices in FigO 

It is now easy to see what role can be played by the 
relevant perturbations in ([3]). These flnite-range in- 
teractions do not need to be very strong to introduce suf- 
flcient mixing between competitive vortex states near the 
pairing transition, and melt the vortex lattice of the su- 
perfluid phase. The obtained normal states are strongly 
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FIG. 3: Density plot of AjT — .Ttl — ^Ah, where JT are 
free energy density minimums of the Abrikosov lattice (AL) 
and a flux lattice (FL) with rtx = 1, Uy = 2, 5q = \/B 
(v / (Amuoa^) — 0.6, T — 0, N = oo, the lowest bosonic Lan- 
dau level only). A.7-" > inside the paired regions bounded 
by the contours, but the halos around the bright domes have 
AT/{2Nlj) < 0.0005. 

correlated vortex liquids. Their properties are not uni- 
versal in the sense that the interactions w„ that stabilize 
them reflect microscopic details of inter-atomic poten- 
tial. Similarly, these perturbations control the properties 
of vortex cores in the superfluid. Unitarity looses many 
universal aspects, but paves the way for the emergence 
of stable vortex liquid states. 
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